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Abstract. We study the production of same-sign W boson pairs at the LHC in double parton interactions. 
Compared with simple factorised double parton distributions (dPDFs), we show that the recently developed 
dPDFs, GS09, lead to non-trivial kinematic correlations between the W bosons. A numerical study of the 
prospects for observing this process using same-sign dilepton signatures, including W ± W ± jj, di-boson 
and heavy flavour backgrounds, at 14 TeV centre-of-mass energy is then performed. It is shown that a 
small excess of same-sign dilepton events from double parton scattering over a background dominated by 
single scattering W ± Z( r y*) production could be observed at the LHC. 



1 Introduction 

The Large Hadron Collider (LHC) at CERN will offer 
many interesting tests of Standard Model (SM) physics 
and, indeed, a precise knowledge of SM processes is neces- 
sary to identify unambiguously New Physics (NP) signals 
from SM backgrounds. The high LHC collision energy will 
open up a new kinematic regime in which certain SM pro- 
cesses will become precisely measureable for the first time. 
An example is multiple parton hard-scattering, i.e. events 
in which two or more distinct hard parton interactions 
occur simultaneously. The theoretical study of such pro- 
cesses goes back to the early days of the parton model [TH3] 
with subsequent extension to perturbative QCD [fflTS], 
Experimental evidence for double parton scattering (DPS) 
has been found in y/s — 63 GeV pp collisions by the AFS 
collaboration at the CERN ISR [T7] , and more recently in 
y/s = 1.8 TeV pp collisions by the CDF collaboration [18] 
and y/s = 1.96 TeV pp collisions by the DO collabora- 
tion [E5] at the Fermilab Tevatron. 

In the standard framework for calculating inclusive 
hard-scattering cross sections in hadron-hadron collisions, 
it is assumed that only one hard interaction occurs per col- 
lision (plus multiple soft interactions). This assumption is 
typically justified on the grounds that the probability of a 
hard parton-parton interaction in a collision is very small. 
Thus the probability of having two or more hard interac- 
tions in a collision is highly suppressed with respect to the 
single interaction probability. 

However, as the collider centre-of-mass energy becomes 
larger, we may expect multiple hard parton collisions to 
become more important. This can be understood as fol- 
lows. Consider a final state consisting of the products of 
two hard collisions A and B, where for example A, B — 



W, Z,jj,tt, ... etc. It is commonly assumed that double 
parton scattering cross sections ^ D ab) can ^ e ^proxi- 
mately factorised into the product of two single scattering 
cross sections: 



DPS 
T (A,B) 



m a (A) CF {B) 
2 Ccff 



(1) 



where the quantity m is a symmetry factor that equals 1 
if A = B and 2 otherwise. The factor a e g in the denom- 
inator has the dimensions of a cross section. The reason 
for this is that given that one hard scattering occurs, the 
probability of the other hard scattering is proportional to 
the flux of accompanying partons; these are confined to 
the colliding protons, and therefore their flux should be 
inversely proportional to the area (cross section) of a pro- 
ton. Of course the (A, B) final state can also be produced 
by a single parton scattering, with cross section <?f AB y 
If the masses of the final states A and B are fixed, then 
the collider energy (y/s) dependence of the hadronic cross 
sections is controlled by the x dependence of the parton 
distribution functions (PDFs) . Because the PDFs increase 
rapidly with x as x — ► 0, the cross sections increase with 
y/s. This implies that crP A B ^ , which is proportional to a 
product of single scattering cross sections, will increase 
more rapidly with y/s than cr^ AB y This raises the possi- 
bility that multiple hard parton scattering, which has re- 
ceived relatively little attention to date, could provide im- 
portant backgrounds to NP signals at the LHC [20-23] . To 
take a simple example for SM Higgs production, the 'stan- 
dard' irreducible background to associated Z + H (— > bb) 
production is qq, gg — > Zbb. However, there is an addi- 
tional DPS background with A = Z and B = bb. 
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Of course the magnitude of the DPS cross section is 
directly dependent on the size of a e g in Eq. ([lj. The re- 
cent CDF and DO measurements [T51[]j5] at the Tevatron, 
utilising the 7 + 3jet final state with A corresponding 
to 7+jet production and B to dijet production, suggest 
cr ff ~ 15 mb, which is roughly 20% of the total (elastic 
+ inelastic) pp cross section at the Tevatron collider en- 
ergy. The non-perturbative physics that determines <7 c ff 
is not well enough understood to be able to predict its 
value at the LHC. If it is proportional to the total in- 
elastic cross section, one might expect a slightly higher 
value at LHC energies. Given this uncertainty, the best 
approach is clearly to find a benchmark double scattering 
process for the LHC, from which a e g can be determined. 
This will not only serve to calibrate DPS backgrounds to 
NP processes, but will also provide important information 
on the non-perturbative structure of the proton. 

In fact, Eq. (TTJ) is only an approximation to the follow- 
ing more general expression: 



DPS 
a {A,B) 



IF ^ 



dx\dx2 dx'i dx 2 d 2 b 



X (xi , £ 2 , b; t\, t 2 ) r k i (x\ ,x' 2 ,b;tx,h) 
xaf k (x 1 ,x' 1 )af l (x 2 ,x 2 ). (2) 

The rij(x\,X2,b;ti,t2) represent generalised double par- 
ton distributions. They may be loosely interpreted as the 
inclusive probability distributions to find a parton i with 
longitudinal momentum fraction x\ at scale t\ = \n(Q 2 ) 
in the proton, in addition to a parton j with longitudinal 
momentum fraction X2 at scale t 2 = ^(Q 2 ), with the two 
partons separated by a transverse distance b. The scale t\ 
is given by the characteristic scale of subprocess A, whilst 
t 2 is equal to the characteristic scale of subprocess B. 

It is typically taken that x 2 , 6; ti, £2) may be 

decomposed in terms of longitudinal and transverse com- 
ponents as follows: 

r ij (xi,X2,b;t 1 ,t 2 ) = D^(xi,x 2 ;ti,t 2 )F-(b). (3) 

Making the further assumption that Fj(b) is the same 
for all parton pairs ij involved in the DPS of interest, this 
leads to: 



OPS 
\A,B) 



2(7, 



cir 



E 

i,j,k,l 



dx\dx 2 dx x dx 2 



x D$ (x 1 ,x 2l t 1 ,t 2 ) Df (x[ ,x' 2 ;t 1 ,t 2 ) 
xaf k (x 1 ,x' 1 )af l (x 2 ,x' 2 ), 



CcfF = 



d 2 b{F{b)f 



(4) 



If one ignores longitudinal momentum correlations such 
that the components DV take the form DV {x\,x 2 ',t\,t2) 
= D\(xx\ ti)D J h (x 2 ] t 2 ), then one finally arrives at the form 
of Eq. . This is the approach that has been taken in ex- 
isting phenomenological calculations regarding DPS [2"0T - 
I27j . Such an approximation is typically justified at low x 



values on the grounds that the population of partons is 
large at these values. 

On the other hand, a number of theoretical studies 
[2"5H5T] have suggested that non-negligible longitudinal 
momentum correlations do exist in the double parton dis- 
tributions (dPDFs) (xi,x 2 ;ti,t 2 ). These papers have 
investigated the special case in which the two factorisation 
scale arguments of the dPDFs are set equal t\ = t 2 = t. 
In [IMJJMI], it is shown that pQCD ('double DGLAP') 
evolution causes the dPDFs to deviate from factorised 
forms, such that even if factorised forms are a good ap- 
proximation at low scales, they cannot be so at higher 
scales. Ref. goes further, and shows that as a result of 
the sum rules the dPDFs have to obey, the PDF factori- 
sation hypothesis D tJ = D l D : > cannot hold for any i,j at 
any scale t - although it may be a reasonable approxima- 
tion for sea quark and gluon distributions at small x. 

Thus, a more accurate way to model DPS is to use 
Eq. (0} along with a set of dPDFs which incorporates the 
effects of pQCD evolution and sum rule constraints (e.g. 
GS09 [29]). The practical implication of the longitudinal 
momentum correlations inherent in the GS09 set is that 
the final states A and B will necessarily be correlated 
in longitudinal momentum, and in particular the rapidity 
distribution of A will not be independent of that of B. 

In this study we focus on same-sign W pair produc- 
tion at the LHC as the paradigm DPS benchmark pro- 
cess, i.e. A = B = W + and A = B = W~ . In this con- 
text, W ± W ± production was first discussed in Ref. [23] 
and subsequently studied in more detail in Ref. [21] and 
Rcfs. [25,26 , the latter of which also compared the cross 
sections of DPS W^W^ and W^W^jj in single scatter- 
ing with the inclusion of kinematic cuts. We will investi- 
gate not only the magnitude of the DPS cross sections, 
but also the rapidity correlations for the DPS W^W^ 1 
final state, both for the GS09 dPDFs and for a simple 
factorised PDF model. Note that the recent CDF and 
DO measurements are not accurate enough to distinguish 
these: the CDF result |T8 j in particular shows no sign of 
x-dependence in the a„Q measured. 

In assessing suitable benchmark processes for DPS at 
LHC, one should of course choose a signal channel for 
which the single scattering background is suppressed. We 
will be interested in W^W^ final states in which both W 
bosons decay leptonically, W — ?> Iv with I = e, /x. This re- 
sults in same-sign dilepton (SSL) signals which we will in- 
vestigate in detail. M /± W ± production has the advantage 
that same-sign single scattering is forbidden at the same 
order in the SM, i.e. there is no ij —> W + W + ,W~W~ 
contribution (cf. qq —> W + W~ , ZZ). The lowest order 
'background' process is W^W^jj, which is of order 0(a 4 ) 
or 0{a 2 s a 2 ). As we shall see, the presence of extra jets 
serves as an efficient tag to veto this background. The 
possibility of using lepton pseudorapidity distribution to 
enhance the DPS signal was considered in [55] [33]. For 
these reasons, it was believed that the same-sign dilep- 
ton channel provides a 'clean environment' for studying 
DPS processes. However there are other important back- 
grounds, for example W^Z{^*) and bb production, both 



Jonathan R. Gaunt et al.: Same-sign W pair production as a probe of double parton scattering at the LHC 



3 



10' 



10' 



i 1 r 



10 2 


"w + ■; 


: : : " ■' " 




: ° w ~ 










10° 


Tevatron 


LHC - 








10" 2 




















,DPS 

w+w- 

T DPS 
W + W + 

w-w- 







4 6 8 10 12 14 16 
CM energy (TeV) 



Fig. 1. Cross sections of various processes in proton (anti-)proton collisions as a function of yfs. The dotted curves correspond 
to single scattering processes, while the solid curves correspond to double scattering processes computed using GS09 dPDFs. 
From top down, the pp single scattering cross sections in the Tevatron region are o w ± , <Jw+w- t (J w ± zj °zz and <y\v ± w ± jj- 
The double scattering cross sections are cr^,+'^ / _ and o~^± s w ± . In the LHC region, the dotted curves are pp cross sections of 



a W+ i a w- > °"w+w-i a W+Z> °zz, &w-z> 



<?w+w+a and a w-w- jj, while the solid curves are (Jw+w- 



DPS a „j „DPS 
'W+W+ U "w-w 



of which can produce a pair of same-sign leptons and miss- 
ing Et- The former leads to a SSL signal when the 'wrong' 
sign lepton from Z decay falls outside the detector accep- 
tance, whereas for the latter, SSL events result when a 
neutral B° meson undergoes B Q -B° mixing, followed by 
leptonic decays. These backgrounds have not been studied 
in detail in the context of di-boson production (a brief dis- 
cussion of the W ± Z(7*) background can be found in [34] h 
We will go beyond comparing the DPS W^W^ signal and 
SPS W^W^jj background p21H2rj] to explore the impact 
of a fairly standard choice of lepton cuts on SSL events 
from both the signal and backgrounds. 

This paper is organised as follows. In the following 
section we discuss in detail the calculation of the sig- 
nal W^W^ — » V^V^vu DPS process cross section. We 
compare total cross sections and rapidity distributions 



obtained using GS09 dPDFs and approximate factorised 
dPDFs. In Sections [3] and S] we study a number of impor- 
tant background contributions to the SSL signal, and in- 
vestigate to what extent they can be suppressed through 
final-state cuts. Our conclusions on the observability of 
DPS at the LHC in the WW SSL channel are presented 
in Section [5j 



2 Signal processes: leptonic channels of 

W ± W ± 

The DPS signal consists of two same-sign leptons and 
missing energy, coming from the decay of two same-sign 
W bosons. The leptons are produced in two simultaneous 
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partonic processes 

ij -> W ± + X l± + vi + X. 



(5) 



The predictions for the DPS signal are calculated accord- 
ing to Eq. (01 . Total cross sections and distributions are 
obtained for four sets of dPDFs: GS09 from [25], and fac- 
torized dPDFs of the fornQ 

Df( Xl ,x 2 ,t) = Dl(x u t)D b h (x 2 ,t)6{l- Xl -x 2 ) 



x(i — xi - x 2 y 



" = 0,1,2, (6) 



hereafter referred to as MSTW n sets. The sPDFs D{{x, t) 
are taken from the MSTW 2008 LO set. As pointed out 
in [25], these factorised sets do not satisfy dDGLAP evo- 
lution or consistent sets of sum rules, although with an 
appropriately chosen value of n they can provide a rea- 
sonable approximation to the 'exact' GS09 dPDFs. The 
factorization scale in the calculations of the DPS signal is 
fixed at \if = Mw for all parton sets. 

The partonic cross sections in Eq. ^ are calculated 
at leading order. At this level, the transverse momen- 
tum (jpt) distribution of the produced W bosons is zero. 
When we come to consider signals and backgrounds with 
realistic experimental cuts, we will need to provide real- 
istic descriptions of kinematic distributions, in particular 
of transverse quantities like the leptonic pr and see 
Section g] 

All numerical results are evaluated with the following 
electroweak input parameters [36]: Mz = 91.188 GeV, 
M w = 80.398 GeV, G F = 0.116637 x 10~ 5 GeV" 2 , T z = 
2.50 GeV and T w = 2.14 GeV. Other EW couplings are 
derived using tree level relations. The effective branching 
ratio BTZ(W + — > using JW as the total width, 

is 0.106. The CKM mixing parameters used are \V US \ = 
0.226, \V ub \ = 0.004, \Vcd\ = 0.230, \V cb \ = 0.041, and 
other parameters are obtained using unitarity constraints. 
We take the value of — 14.5 mb, consistent with the 
Tevatron measurements. 

A comparison of the cross sections for DPS and SPS 
WW production processes as a function of collider centre- 
of-mass (CM) energy y/s is shown in Fig. [1] The DPS 
cross sections are obtained with the GS09 set of dPDFs. 
We see immediately that while the single scattering qq — > 
W + W~ , W ± Z 1 ZZ cross sections dominate at all collider 
energies, the DPS W ± W ± and SPS W ± W ± jj cross sec- 
tions are comparable in magnitude. As we shall demon- 
strate, with an appropriate set of jet veto cuts the SPS 
W^W^jj cross sections can be significantly reduced. We 
quote the values of the total cross section for DPS WW 
production processes for various CM energies at the LHC 
in Table [1] In Tables Q] and [2J we show also the values of 
the ratio „ 

a w+w- 

which measures the deviation from the factorisation ap- 
proach, Eq. ([1]). When R — 1, factorisation is exact. We 
note that factorisation is broken at the 20% to 30% level, 
and the approximation improves at the higher collider en- 
ergies as lower x regions are probed. 





O"GS09 




s/s =7 TeV 


yfs = 10 TeV 


s/s = 14 TeV 


w+w~ 


0.107 


0.250 


0.546 


W+W+ 


0.0640 


0.148 


0.321 


w~w~ 


0.0317 


0.0793 


0.182 




R 




0.709 


0.751 


0.784 



Table 1. DPS WW total cross sections (in pb) for pp colli- 
sions at different CM energies. The values are obtained by first 
calculating the leptonic cross sections, before dividing by the 
corresponding branching ratios. All cross sections are evalu- 
ated using GS09 parton distributions. The ratio R measures 
deviation from the factorisation approach, as explained in the 
text. 





"GSOO 


"MSTWn 


"MSTW, 


CT MSTW, 


w + w~ 


0.546 


0.496 


0.409 


0.348 


W+W+ 


0.321 


0.338 


0.269 


0.223 


w~w~ 


0.182 


0.182 


0.156 


0.136 




R 




0.784 


1.00 


1.00 


1.00 



Table 2. DPS WW total cross sections (in pb) for pp collisions 
at y/s = 14 TeV evaluated using different dPDFs sets. The val- 
ues are obtained by first calculating the leptonic cross sections, 
before dividing by the corresponding branching ratios. 



Since the DPS total cross sections are most relevant 
for the LHC with y/s — 14 TeV, we will focus our anal- 
ysis on this case. It is interesting to investigate the effect 
of the different sets of the dPDFs, cf. Table [2] Note that 
all the factorised distributions contain a kinematic factor 
9(1 — X\ — x 2 ), which automatically implies R ^ 1 for these 
sets, although the effect is numerically very small at LHC 
(14 TeV) energies where small x values are probed. The 
DPS total cross sections for same-sign W ± W ± production 
are very similar when using the GS09 and MSTWo sets 
whereas the values for the MSTWi and MSTW 2 sets are 
smaller due to the suppression caused by the (1 — x\ — x 2 ) n 
factors. The DPS opposite-sign W ± W T cross section is 
somewhat larger for GS09 than for MSTWo- This is be- 
cause of an enhancement of the D q ^ (q — u.d) GS09 
dPDFs induced by the g — > qq contribution to dDGLAP 
evolution [29j. We also observe that the MSTW W + W + 
cross section is slightly larger than its GS09 counterpart, 
whereas such a difference is not observed in the W~W~ 
cross sections. The probable explanation behind this is 
the following. Both the uu and dd inputs to GS09 consist 
of (roughly) factorised forms minus terms to take account 
of number effects (^D Uv D Uv and D d "D d " respectively). In 
absolute terms, the subtraction from the uu input is larger 
than that from the dd (since ±D u *D Uv 2D d "D d -). This 
results in a greater decrease of the GS09 W + W + cross sec- 
tion relative to that of MSTWo than occurs for W~W~ 
production. In fact, the W~W~ cross section for GS09 
ends up being roughly equal to that of MSTWo, since all 
GS09 dPDFs are 'fed' by the sPDFs during evolution (ei- 



We use the notation t = In fip . 
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Fig. 2. Normalised lepton pseudorapidity distributions for pp collisions at y/s = 14 TeV evaluated using different dPDFs. No 
cuts are applied. 
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Fig. 3. Pseudorapidity asymmetry a m for pp collisions at \/s = 14 TeV evaluated using different dPDFs. No cuts are applied. 



ther directly or indirectly via the gg dPDF). The 'sPDF 
feed' contribution to the GS09 dPDFs compensates for the 
subtraction from the dd input in the case of the W~W~ 
cross section. On the other hand, a similar sPDF feed con- 
tribution cannot counterbalance the greater subtraction 
from the uu input in the W + W + case, causing the GS09 
W + W + cross section to be smaller than the MSTWo cross 
section. Obviously, the deviations from the factorization 
approach measured by the ratio R are small when using 
the MSTW„ sets. In GS09, there are sizeable deviations 
from factorisation, particularly in the valence sector (i.e. 
jjgvqv D qv D qv ), and these are reflected in smaller val- 
ues of R. 

In order to study the potential of using same-sign WW 
production for measuring DPS, a full study of the sig- 
nal, including leptonic decays and cuts, is necessary. In 
Fig. 2(a) wc present the normalised pseudorapidity distri- 
bution of the l + leptons, for different sets of dPDFs. 



Although visible, the differences in the lepton pseudo- 
rapidity distributions for the various dPDFs are small. 
The effect of the correlations in the longitudinal momen- 
tum fractions in the GS09 set is closely reproduced by 
the MSTWi set, a result which could be expected from 
the compar ative s tudy of GS09 and MSTW„ distributions 
in [33]. Fig. 2(b) shows the corresponding plot for l~ lep- 
tons. Although the shapes of the pseudorapidity distri- 
butions are different (a simple reflection of the difference 
between u and d quark PDFs) , the qualitative differences 
between the various sets for l~ are the same as for l + 
production. 

The sensitivity to longitudinal correlations can be max- 
imized in the following asymmetry 



ojjTh x Vh < 0) ~ <?(Vh x r] h > 0) 
v(Vh x Vh < 0) + cr(r/ h x r] h > 0) ' 



(7) 
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where r\\ is the lepton pseudorapidity, and ^(J, |t?z 2 | > 
?7 ; mm . The asymmetry measures the extent to which the 
presence of one W produced at high rapidity affects the 
probability of finding another W boson with similarly 
large rapidity. A positive a Vl means that the leptons pre- 
fer to lie in opposite hemisphere. For higher values of ?7 ; mm 
the effec t of t he correl ations becomes more pronounced, 
cf. Figs. 3(a) and |3(b)| Such behaviour is to be expected 
as the correlations are most important for the distribu- 
tions probed at high values of x for both partons in the 
same proton, reached when the leptons (W bosons) are 
produced at high rapidities. 

In this section we have restricted our attention to basic 
quantities (total cross sections and pseudorapidity distri- 
butions) to illustrate the impact of the various dPDF sets, 
and in particular the effect of the correlations that are a 
feature of the GS09 set. As we shall see below, in practice 
we have to introduce further cuts on the final-state parti- 
cles in order to suppress large backgrounds from SPS pro- 
cesses. In Section @] we will investigate the effects on DPS 
event rates and correlations after imposing these cuts. 
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Fig. 4. a(pp -> W + W + jj) 
function of max prj and T?™ ln 



[BTZ(W+ -> fi+^)} 2 (fb) 
No other cuts are applied. 



3 Backgrounds 

We now turn to discuss background contributions. Re- 
call that a signal event consists of two same-sign charged 
leptons plus missing transverse energy, where the leptons 
originate in two separate hard scatterings from one pp col- 
lision. As already mentioned, W+W+jj single scattering 
provides the lowest order, irreducible background. Heavy 
flavour production and gauge boson pair production can 
also lead to same-sign dilepton events. We shall discuss 
the relevant properties of these background contributions, 
and propose cuts that can reduce them. 



example, we see that vetoing events with jet pr > 20 GeV 
with no pseudorapidity requirement suppresses the cross 
section to about two orders of magnitude below the signal 
(DPS) cross section. We conclude that this background 
can be effectively suppressed by a jet veto, and so will not 
be discussed further. However we note in passing that this 
process may contribute to SSL backgrounds in beyond SM 
scenarios, particularly in event topologies including jets, 
for example single slepton production in supersymmetry 
without R-parity. 



3.2 The heavy flavour background 



3.1 The W+W+jj single scattering background 

In single scattering, the lowest order in which same-sign 
W's can be produced are C(a|a 2 ) and 0(a 4 ). By consid- 
ering overall charge conservation, one sees that all the par- 
ticipating partons must be quarks. The QCD diagrams in- 
volve two quark propagators and exchange of a t— channel 
gluon, whereas neutral electroweak gauge boson and Higgs 
exchanges are also possible for the EW contributions. As 
shown in Fig. [TJ the cross sections for W+W+jj and the 
DPS W+W+ signal are of the same order of magnitude 
at LHC energies. 

The presence of jets provide an important handle to 
suppress this single scattering background. The presence 
of two W's implies that the jet energies are typically of 
the order of the weak scale, Myy, and their px distribu- 
tion can have a long tail. A veto of events with (central) 
high pt jets will therefore be useful in suppressing this 
background. In Fig. 01 we show the variation of the cross 
section of <r(pp -> W+W+jj) ■ [BK(W+ ^+v^)] 2 with 
maximum transverse momentum (max pxj ) and minimum 

pseudorapidity (rff 1111 ) of the jets allowed in an event. For 



The heavy flavour production processes pp —> QQ + X 
with Q — t,b can lead to SSL events. In ti production, 
the dominant contribution to same-sign dilepton events is 
where both a top and the bottom of the other top decay 
semi-leptonically, for example 



t-> W+b^l+vb, 
i -> W~b-> qq'l+pc. 



(8) 



The final state is therefore nominally 1+1+ + 2v + 4 jets. 
Requiring two leptons of the same-sign means that one of 
the leptons should come from b— quark decay. This is im- 
portant because the b— quark, originating from top decay, 
is energetic and its daughter lepton will generally not be 
isolated from the other (hadronic) decay products. Hence 
a tight lepton isolation requirement will be effective in re- 
ducing this background. The lepton from the W typically 
acquires a harder px spectrum than the signal leptons. 
Imposing a maximum lepton pt cut will therefore also be 
useful. Similarly to the W+W+jj process, the presence 
of energetic jets again provides an effective suppression 
mechanism. Because of these considerations, we expect 
that ti production will only contribute subdominantly to 
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the background, and therefore will also be neglected in our 
numerical analysis. 

In bb production, the hard process is initiated by two 
gluons. To obtain same-sign dileptons, the two resulting 
B mesons should decay semi-leptonically. Further, one of 
the B mesons must undergo B°-B° mixing before decay, 
which is possible if the B meson is neutral. The relevant 
processes are thus: 

gg -> bb^r BB + ... , 
B l + vX, 

B° -> B° -> 1+vX, (9) 

together with the charge conjugation processes. 

Due to its large production cross section at the LHC, 
this can be a problematic background (a detailed study 
of heavy flavour backgrounds to lepton-pair-plus-missing- 
transverse-energy final states can be found, for example, 
in Ref. |35|). Fortunately, the kinematic properties of this 
process are very different from those of the signal, as the 
two scales involved are very different (mt vs. Mw)- For 
example, the lepton pt from bb decay peaks at very low 
values and decreases exponentially as pt increases. As the 
Pt of the leptons comes primarily from the pt of the par- 
ent B mesons, at high pt the lepton will tend to align with 
other hadronic decay products making isolation difficult. 
An energetic isolated lepton will also tend to be accompa- 
nied by a soft neutrino, leading to low transverse missing 
energy $t- As a result, a combination of lepton isolation, 
minimum px and $t cuts should effectively suppress this 
background. We shall study the impact of these cuts in 
the numerical simulation in the following section. 

The cc process also leads to same-sign dilepton events 
through a mechanism similar to that of bb. In this case, the 
lepton pt spectra peaks at even lower values. As well as 
having a strongly suppressed contribution in our px region 
of interest, the leptons are again difficult to be isolated 
from the other decay products, and low $t is expected. 
As a result, the lepton cuts used to suppress the bb back- 
ground should effectively suppress cc events as well. 

3.3 The electroweak gauge boson pair background 

In principle, the production and leptonic decay of heavy 
weak boson pairs, V V —> 4 leptons with V = W,Z can also 
provide a sizeable background. However these processes 
do not naturally give rise to events with same-sign lepton 
pairs as the only visible particles. For example, W Z(y*) 
and Z(7*)Z(7*) can lead to 3 or 4 charged leptons re- 
spectively, and can mimic the same-sign di-lepton signal 
if the 'wrong' sign leptons are not detected. This happens 
if they fall outside the detector acceptance, or if they are 
not reconstructed. The relevant processes are then 

qq' -> W+Z(-f*) -> l+vl+(r), 

qq ->Z( 7 *)Z( 7 *) -►/+(/-)/+(*-), (10) 

and their charge-conjugated processes. In the above ex- 
pressions the leptons in brackets are not identified. Clearly, 



a wrong sign lepton veto will be able to reduce this back- 
ground when more than 2 leptons are identified. The vir- 
tual 7* can still contribute significantly when it decays 
asymmetrically into a hard and a soft lepton in the cen- 
tral region. Following [37| , this may be suppressed by look- 
ing for isolated charged tracks that form a low invariant 
mass with one of the same-sign leptons. For the Z, the 
dominant contribution is when the wrong sign lepton lies 
outside the central region and is not reconstructed as a 
lepton. This in turn pulls its partner lepton towards the 
large pseudorapidity region, potentially providing a shape 
variable for further discrimination. Furthermore, the lep- 
ton px spectrum extends beyond that of the DPS signal, 
and therefore a maximum lepton pt cut will be useful. 

3.4 Other backgrounds 

Another possible source of same-sign lepton pairs comes 
from multi-particle interactions, when production of W's 
of the same sign from two separate proton collisions occurs 
during the same bunch crossing. In the factorised dPDF 
approximation, this background and the signal, where the 
double scattering occurs in one proton-proton collision, 
are expected to exhibit similar kinematic properties. How- 
ever the leptons from the multi-particle-interaction back- 
ground will in general have tracks pointing back to two 
different locations along the beam axis. In the Appendix, 
we perform a simple estimation to show that reasonable 
longitudinal vertex resolution in the LHC detectors should 
allow significant suppression of these events. 

Apart from the physics backgrounds discussed, non- 
physics backgrounds can also be important. For example, 
W+ jets may contribute when a jet is mis-identified as a 
lepton, or ZZ production may contribute when one Z de- 
cays invisibly and the charge of a lepton from the decay of 
the other Z is mis-identified. A thorough investigation of 
these effects requires a detailed detector simulation. This 
is beyond the scope of the present study, and so we restrict 
ourselves to studying the effect of physics backgrounds 
only. 

4 Numerical Study 

In the previous two sections we have discussed the DPS 
signal, and its characteristic properties, and the most im- 
portant single scattering backgrounds. In this section we 
carry out numerical studies to investigate the relative size 
of the various contributions and how to reduce them. 

We first discuss the cuts, referred to as 'basic cuts' be- 
low, applied in our numerical analysis. They are 'basic' 
in the sense that they are necessary to reduce the back- 
ground to a manageable level, while keeping the signal 
largely intact. Later we shall discuss refinements that can 
improve the signal to background ratio, but at a price of 
further reducing the signal. 

As discussed in the previous section, the W ± W ± jj 
and tt will (the latter partly) be suppressed by a jet veto. 
In particular, rejecting events with jets having pt > 20 
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GeV will suppress the W ± W ± jj background to about two 
orders of magnitude smaller than the DPS signal before 
application of further cuts. A jet veto is thus implicitly 
assumed in the following, and only the bb, W ± Z(~f*) and 
Z(-f*)Z(-f*) backgrounds are simulated. 
The basic cuts are as follows: 

1. Both lcptons in the like sign lepton pair must have 
pseudorapidity \r]\ < 2.5. 

2. Both leptons are required to be isolated: 

-E^so — ^So = 10 GeV, where E l ISO is the hadronic 
transverse energy in a cone of R = 0.4 surrounding 
each of the like-sign leptons. 

3. The transverse momenta of both leptons, p l T , must sat- 
isfy 20 < p l T < 60 GeV. 

4. An event is rejected whenever a third, opposite-signed, 
lepton is identified. A lepton is assumed to be identified 
with 100% efficiency when p l T > pl^ and \rj\ < i] id , 
where pif, = 10 GeV and i] id = 2.5. 

5. The missing transverse energy Tfix of an event must 
satisfy # T > 20 GeV. 

6. Reject an event if a charged (lepton) track with pi£ > 
px > 1 GeV forms an invariant mass < 1 GeV with 
one of the same-sign leptons. 

In the above, the muons are treated as invisible particles 
when 1 77 1 > 2.5. We define Jpx to be the magnitude of the 
vector px sum of all visible particles. Note that the above 
cuts are designed to be within the capabilities of the LHC 
general purpose detectors, ATLAS and CMS. All our cross 
sections are evaluated at y/s = 14 TeV. 

We next discuss event simulations. The bb events are 
generated using HERWIG6.510 [38] ■ The large total cross 
section a bb ~ 500 [A> means that, in practice, a parton- 
level cut on the transverse momentum of the bottom quarks, 
p T , is imposed to make the simulation manageable. This 
parton-level cut is chosen to be p T > 20 GeV, which is 
motivated by one of the basic cuts, p l T > 20 GeV, dis- 
cussed above. The resulting lepton px distribution should 
remain unchanged, as most isolated leptons from B meson 
decays have p l T smaller than the transverse momentum of 
the parent b quark. To improve efficiency of event genera- 
tion, we further force the B mesons to always decay semi- 
leptonically. Whenever one or more neutral B mesons are 
produced, exactly one of them undergoes B°-B° mixing. 
As discussed in [35], this neglects the production of lep- 
tons in charm decays. However these leptons are expected 
to have a lower px and be less well isolated than the lep- 
tons produced in b decays. The amount of Monte Carlo 
data generated is equivalent to about 200 fb . The cross 
section a bb (p T > 20 GeV) is then normalised to the cor- 
responding value obtained in MCFM [40 using the central 
MSTW LO sPDF set. This value is found to be a bb (p T > 
20 GeV) = 5.15 /ib, where the renormalisation and fac- 
torisation scales are set at fx R = fj, F = m b = 4.75 GeV 
respectively. 

Simulations of the other processes are performed at 
leading-order parton level only. Given that the value of 
a e ff, which controls the overall magnitude of the DPS sig- 
nal, is unknown at LHC energies, this approximation is 







(fb) 


W^W* (DPS) 


0.82 


0.46 




5.1 


3.6 


bb (p b T > 20 GeV) 


0.84 


0.67 


0.43 


0.43 



Table 3. Cross sections (in fb) of the processes simulated af- 
ter the basic cuts are applied, including branching ratios corre- 
sponding to same-sign dimuon production. The cross sections 
for n + pb + and yT yT production are displayed in the second 
and third columns respectively. The signal cross sections are 
obtained assuming <r e g- = 14.5 mb. 



sufficient for the present purpose. The parton distribu- 
tions and electroweak parameters used are discussed in 
Section [21 The value of as(Mz) is 0.13939, and we use 
one loop as running throughout. 

For W ± Z(j*) and Z(j*)Z(j*), we obtain the matrix 
elements from MADGRAPH gTJ[22], including both doubly 
and singly resonating diagrams. The phase space integra- 
tion is performed using VEGAS [43] , All gauge bosons are 
decayed leptonically, and off-shell effects are included in 
the simulations. The factorisation scale is chosen to be 
[if = Mw for all these (leading-order) processes. 

The signal processes W ± VK ± are generated in a similar 
way. Here GS09 is used, and again we set \ip = Mw- The 
effective cross section cr e g- is taken to be 14.5 mb, its value 
obtained by the CDF collaboration at the Tevatron. Since 
we use cuts on the px of the leptons, it is important to 
generate realistic px distributions. To account for the non- 
zero px of a W boson, we introduce a 'p^-smearing'. It is 
well known that the fixed-order perturbative calculation 
fails to describe the px distribution of single gauge bosons 
produced at small px and that the correct description 
is provided by the resummed calculations, supplemented 
by a parameterization of non-perturbative effects at very 
small px- Resummation takes into account modifications 
to the px distribution due to multiple soft gluon emission. 
We calculate the resummed distribution to next-to-leading 
logarithmic (NLL) accuracy with NLO MSTW pdfs, us- 
ing the code of Ref. |44j and the non-perturbative param- 
eterisation of Ref. 05] ■ The leptons originating from the 
two (virtual) W bosons are then boosted independently 
in the azimuthal plane according to the px distribution 
obtained in this way. Even though the dPDFs used in this 
study are leading-order quantities, which means that for- 
mally it is more consistent to adopt a leading logarithmic 
(LL) calculation, we argue that the described (NLL) pro- 
cedure is more appropriate for our studies, since it gives a 
more realistic description of the W px distribution. How- 
ever it should be kept in mind that different 'px smearing' 
prescriptions will lead to slight changes in kinematic dis- 
tributions. 

The cross sections for the simulated DPS signal and 
background processes are displayed in Table [3J The cross 
sections correspond to detecting same-sign dimuon events 
only. Inclusion of e ± e ± and e ± y ± events can be estimated 
by multiplying the results in Table [3] by appropriate fac- 
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Fig. 5. Selected kinematic distributions for the positively charged (H — h) SSL events. The DPS signal, W + Z(j*), Z('y*)Z{'y* 
and bb backgrounds are described by solid black, green dotted, blue dotted and red dashed lines respectively. 



torsQ We see that although the total bb cross section is 
huge, the basic cuts are very effective in suppressing this 



2 Due to different dPDFs, cuts and lepton channels consid- 
ered, our numerical values for the cross sections are different 
from those presented in Ref. [26] . 



background. Instead the largest background comes from 
WZ(j*) production, whereas the Z("f*)Z(~f*) background 
is about a factor of 6 smaller. In the latter, we observe an 
asymmetry of the ^+ and cr^-^- cross sections. These 
two processes would have the same cross sections if the 
parton level amplitudes were the same upon interchang- 
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ing the momenta of an opposite-sign lepton pair. However 
this not the subset of the Feynman diagrams 

contributing to these processes changes sign, leading to 
different kinematic (in particular pseudorapidity) distri- 
butions 

One could consider whether kinematic distributions 
might be able to further discriminate the signal from back- 
ground. Given that the signal cross sections are only in the 
1 fb region, care is required not to significantly suppress 
the signal. 

In Fig.[5]and Fig.[S]we show some representative kine- 
matic distributions for the positively charged (++) and 

negatively charged ( ) SSL events respectively. In these 

figures, we use only ten bins in each plot to reflect the an- 
ticipated modest number of events expected for the fem- 
tobarn level processes considered, at least in early LHC 
data taking. 

For the signal, the (++) r\i distribution is almost flat, 
increasing slightly towards large \r)i\, whereas the signal 

for the ( ) events shows a small central peak instead. 

Both of these features are characteristic of charged lep- 
tons from single Drell-Yan production. The distribution of 
A(fiu , the angular separation of the SSL pair in the trans- 
verse plane, dips near A(f>u = 7r, as the neutrinos tend to 
be produced back-to-back in this limit, thus reducing the 
missing transverse momentum $t- 

For the W Z(j*) contributions, small peaks appear in 
the Atfra — and tot/; = regions. Here ttitu is the trans- 
verse mass of the SSL lepton pair. These two peaks are 
correlated, and arise primarily from the 7* contributions. 
The Z(j*)Z(-f*) background exhibits similar peaking fea- 
tures for the same reason. Furthermore, the rji distribution 
for this background peaks strongly towards the high re- 
gion. This is because the Z contributes most significantly 
when one of its daughter leptons is produced outside the 
central tracking region, pulling its partner forward as a 
result. 

Overall, we see that while the bb distributions are dis- 
tinct from those of the signal, the dominant W ± Z( , y*) 
background has kinematic distributions fairly similar to 
the signal. A cut on events when the lepton pair are ap- 
proximately back-to-back in the transverse plane, for ex- 
ample A(f>u < 0.77T, could help, but further kinematic cuts 
are unlikely to improve the signal to background ratio sig- 
nificantly. Optimisation of the cuts would in any case re- 
quire a more detailed detector simulation, and is beyond 
the scope of this work. 

On the positive side, it might be advantageous to ex- 
ploit the fact that the value of a r/l is relatively small for 
the DPS signal. This means that the probability of the 
same-sign leptons from the signal process to lie in oppo- 
site (771 x 7] 2 < 0) and same (771 x r) 2 > 0) hemispheres 
are approximately equal. However for the single scatter- 
ing background, the final states prefer to have small Ar\ 

3 It is interesting to note that this happens even for diagrams 
involving only 7* and not Z. The parton level amplitudes of 
these diagrams do not change under charge conjugation. How- 
ever the associated PDF weights change as a result, leading to 
different kinematic distributions at the hadron level. 



in order to minimise the partonic centre-of-mass energy. 
This implies that a m will generally be negative. In Fig. [JJ 
we show a m as a function of rjf" 1 after the basic cuts are 
imposed. Requiring the leptons to satisfy rjn x r]i2 < 
would reduce the signal by about a factor of 2, while re- 
ducing 2/3 of the W ± Z{ 1 *) and almost all the Z(7*)Z( 7 *) 
backgrounds. 

The ratio of positively charged (++) and negatively 

charged ( ) SSL events may also be used. Due to the 

(quasi-)factorised nature of DPS production, the signal 
cross ratio a^w+/ a w~w- * s ex P ec ted to be similar to 
(a w + / <j w - ) 2 . The former ratio can be inferred from Ta- 
ble [2] to be about 1.8 at 14 TeV, whereas the square of 
a w+l a w- at LO is ~ 1.36 2 = 1.85. An important point 
is that this ratio is stable against the application of cuts, 
as can be inferred from the results in Table [3J Because it 
originates in gg scattering, the bb background contributes 

equally to (H — h) and ( ) lepton pair production. From 

Fig. [51 we see that the <J W + z /cr w - z ratio is roughly 1.4, a 
reflection of the u/d parton distribution ratio in the pro- 
ton. For the Z( r y*)Z(j*) process, the corresponding ratio 
is 1.25 after including the basic cuts. 

Even though these ratios are cut dependent and would 
presumably change slightly after properly including higher- 
order effects, the presence of DPS events predicts that the 
ratio (++)/( ) would exceed the expectation from sin- 
gle scattering background events only. 



5 Summary 

We have studied same-sign W pair production as a poten- 
tial signal of double parton scattering at 14 TeV at the 
LHC. We considered both the DPS signal (in the purely 
leptonic decay channel) and a number of backgrounds 
from the usual single parton scattering. 

We first showed that the improved GS09 dPDFs lead 
to different production properties compared to naive fac- 
torisation models. In particular, the use of GS09 dPDFs 
leads to non-trivial rapidity correlations of the final-state 
leptons. Valence number conservation implies that it is un- 
likely to find 2 high- a; valence quarks of the same flavour 
from a single proton, resulting in a positive lepton pseudo- 
rapidity asymmetry a Vl for the DPS processes. Otherwise, 
the lepton pseudorapidity distribution using GS09 can 
be reasonably well approximated by the MSTWi quasi- 
factorized model. 

Our calculations of the DPS signal were undertaken 
with the assumption that a e Q is a constant, and equal 
to the CDF measured value of 14.5 mb. We remind the 
reader that although this assumption is consistent with 
existing experimental data, it is likely that there will be a 
certain degree of process dependence in a e g [15] ■ It is also 
possible that different scale factors might be appropriate 
for contributions to the dPDF coming from different terms 
in the dDGLAP equation [21]. Unless the a e g for equal 
sign WW production at the LHC is very much smaller 
than the CDF value, it is likely that the theoretical aspects 
just described will be difficult to study in the early stages 
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Fig. 6. Selected kinematic distributions for the negatively charged ( ) SSL events. The DPS signal, W Z(y*), Z(y*)Z('j* 

and bb backgrounds are described by solid black, green dotted, blue dotted and red dashed lines respectively. 



of LHC running due to limited statistics. However, the 
fundamental characteristic of the DPS signal process - 
two (quasi-) independent scatterings - means that rjn x r/i2 

and the (++)/( ) event ratio could provide additional 

experimental handles. 



On the other hand, more work needs to be done in 
order to suppress the background further. In addition to 
the W^W^jj background previously considered in the 
literature, we have also considered di-boson and heavy 
flavour backgrounds in some detail. Even though the 'non- 
reducible' W^W^jj single scattering background can be 
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be identified. From bottom to top, the lines correspond to =5, 10 and 15 GeV respectively. In Fig. 8(b) the ratio for the 
DPS signal is also shown. 



effectively suppressed by a central jet veto, the di-boson 
background remains significant after applying a fairly ba- 
sic set of cuts. Using the canonical a e g — 14.5 mb to 
calibrate the DPS signal, we find that the diboson back- 
ground can be a factor of a few larger than the signal. 

Given the exploratory nature of this study, we have 
not included cither detector effects or higher-order per- 
turbative corrections. Given the similarity between the 
W ± Z(7*) and signal distributions, it might not be very 
beneficial to impose simple additional physics cuts. How- 
ever, the Standard Model processes leading to tri-lepton 
events will be studied in detail at the LHC. An in-depth 



understanding of this process might allow accurate extrap- 
olation into our region of interest, where one of the leptons 
is not detected. This would provide an important calibra- 
tion of the diboson background. Furthermore, 'detector 
cuts' could also be made. For example, tighter lepton iso- 
lation and wrong-sign lepton vetos could be used. Extend- 
ing the pseudorapidity region where a (wrong sign) lepton 
can be identified could also be useful. Displaced charged 
lepton vertices might also be helpful in further reducing 
the bb background. 

In summary, our preliminary 'first look' analysis has 
shown that a small excess of SSL events from double par- 
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ton scattering could be observed at the LHC, and we are 
optimistic that further improvements can be made to en- 
hance the signal. 
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Appendix: multiple particle interaction 

Given a luminosity (L), a single scattering cross section 
(a) and rate of bunch crossing (£?), the average number 
of events per bunch crossing, (n) , is given by 



La 
~~B' 



(11) 



Using Poisson's statistics, this can be translated into a 
multiple particle interaction cross section, aN, where 



-<n) 



(n) N B 



Nl L 
~N!\BJ 



(12) 



where the last equation holds if (n) < 1. An effective mul- 
tiple particle interaction parameter, a N e g, can be defined 
analogous to a e g for multi-parton interaction: 



a N 



Veff 



(13) 



The above approximation ceases to be valid when a ~ 
a n eff' wn i cn implies also that (n) ~ 1. For B = 4 ■ 10 7 s _1 



and L = 10 34 cm 2 s x , we get a 



jv,eff 



4 mb. However, 



in double particle interactions, the two interaction vertices 
typically do not overlap. Using the RMS bunch length 
of 7.5 cm |46j and intrinsic z-resolution of 115 /mxi and 
580 /im for the Pixel detector and SCT at ATLAS g7] re- 
spectively, the probability that 2 independent scatterings 
overlap each other is estimated to be of the order O(0. 1)%. 
As a result this background contribution to double parton 
scattering, assuming a e Q = 14.5 mb, can be safely ne- 
glected. 
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